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We investigate a possible nuclear correction to the NuTeV measurement of the weak-mixing 
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In particular, a nuclear modification difference between u„ and d v distributions 



contributes to the NuTeV measurement with the iron target. First, the modification difference is 
determined by a \ 2 analysis so as to reproduce nuclear data on the structure function F2 and Drell- 
Yan processes. Then, taking the NuTeV kinematics into account, we calculate a contribution to the 
sin 2 6w determination. In addition, its uncertainty is estimated by the Hessian method. Although 
the uncertainty becomes comparable to the NuTeV deviation, the effect is not large enough to 
explain the whole NuTeV sin 2 dw anomaly at this stage. However, it is difficult to determine such 
a nuclear modification difference, so that we need further investigations on the difference and its 
effect on the NuTeV anomaly. 



PACS numbers: 13.15.+& 13.60.Hb, 24.85.+p 



I. INTRODUCTION 

Weak-mixing angle sin 2 9w is one of fundamental 
constants, and there are various experimental methods 
for determining it pj. They include experiments on 
atomic parity violation, left-right and forward-backward 
asymmetries in electron-positron annihilation, W and Z 
masses, elastic and inelastic neutrino scattering, parity- 
violating electron scattering, and polarized M0ller scat- 
tering. 

The NuTeV collaboration announced that their 
neutrino-nucleus scattering data indicate a significant 
deviation from the standard model. If the neutrino- 
nucleus scattering data are excluded, the average value 
is sin 2 9 W = 0.2227 ± 0.0004 Q in the on-shell scheme. 
The NuTeV value, sin 2 W = 0.2277 ± 0.0013 (stat) ± 
0.0009 (syst) @, is significantly larger than the aver- 
age of the other data. This difference is called "NuTeV 
anomaly" . Independent experiments are in progress to 
measure the angle by the M0ller scattering 4] and parity- 
violating electron scattering ||. It is also interesting to 
find Q 2 dependence of sin 2 0w from various experiments. 
From these measurements, we expect that the experi- 
mental situation will become clear. However, since it is 
an important fundamental constant, fast clarification of 
the NuTeV deviation is needed. 

There are a number of papers which intended to ex- 
plain the deviation. Although there may be new mech- 
anisms beyond the standard model 0, we should inves- 
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tigate all the possibilities within the present theoretical 
framework. Next-to-leading-order effects are not a source 
of the deviation 0. There are proposals for conven- 
tional explanations 0. First, because the used target 
is the iron nucleus instead of the isoscalar nucleon, it 
is natural to explain it in terms of nuclear corrections 
[H El El However, there is still no clear inter- 
pretation at this stage. Second, the difference between 
strange and antistrange quark distributions should con- 
tribute to the deviation El El El El- However, it 
cannot be uniquely determined at this stage. For exam- 
ple, a CTEQ analysis result [15( differs from a NuTeV 
result E3- Third, there are proposals that the anomaly 
could be explained by isospin violation |17|. which could 
be related to the nuclear correction in Ref. as shown 
in Sec. ITvTfl 

In this paper, we investigate a conventional explana- 
tion in terms of a nuclear correction. We study the details 
of the nuclear effect in Ref. El 1 where it was pointed 
out that the nuclear modification difference between the 
u v and d v distributions contributes to the anomaly. Be- 
cause such a nuclear modification is unknown, possible 
differences are simply estimated by considering baryon- 
number and charge conservations |ll| . It is the purpose 
of this paper to clarify this kind of nuclear effect by a x 2 
analysis with many nuclear data by assigning explicit pa- 
rameters for the difference. We also estimate uncertain- 
ties of the difference by the Hessian method, and they 
are especially important for discussing an effect on the 
sin 2 6w anomaly. 

This paper is organized as follows. In Sec. [H] we 
show a nuclear modification effect on sin 2 6w by using 
the Paschos-Wolfenstein relation. Next, a x 2 analysis 
method is explained in Sec. IHII for determining the mod- 
ification difference between u v and d v . Analysis results 
and an actual effect on the NuTeV sin 2 8\y are discussed 
in Sec. IIVI These studies are summarized in Sec. W\ 



II. CORRECTIONS TO 
PASCHOS-WOLFENSTEIN RELATION 

Neutrino and antineutrino scattering data are used for 
extracting the weak mixing angle by the NuTeV collab- 
oration. In discussing various corrections to the NuTeV 
measurement, it is theoretically useful to find the cor- 
rections to the Paschos-Wolfenstein (PW) relation [Isj : 
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where <Jcc an d <j 'nc are charged-current (CC) and 
neutral-current (NC) cross sections, respectively, in 
neutrino-nucleon deep inelastic scattering. 

The various correction factors to the PW relation are 
discussed in Ref. Because they are important for 

our analysis, an outline is discussed in this section. The 
iron target is used in the NuTeV experiment, so that nu- 
clear corrections should be properly taken into account. 
The correction investigated in this paper is the nuclear 
modification difference between the u v and d v distribu- 
tions. Nuclear modifications of the valence-quark distri- 
butions, u v and d v , are expressed by w Uv and w dv , which 
are defined by 



u*{x,Cf) = w Uv {x,Q\A,Z) 
d A (x,Q 2 ) = w d Jx,Q 2 ,A,Z) 



Zu v (x,Q 2 ) + Nd v {x,Q 2 ) 
A 

Zd v {x,Q 2 )+Nu v (x,Q 2 ) 



Here, u A and d A are the up- and down- valence quark 
distributions, respectively, in a nucleus; u v and d v are 
the distributions in the proton. The variable x is the 
Bjorken scaling variable, Q 2 is defined by the momentum 
transfer q as Q 2 = —q 2 , Z is the atomic number of the 
target nucleus, and A is its mass number. In Ref. |19| . 
the weight functions w Uv and w dv are defined only at the 
fixed Q 2 point, Q 2 = l GeV 2 ; however, they are defined 
at any Q 2 in the present research. As mentioned in Ref. 
[20|. the form of Eq. J2J is not appropriate at x — ► 1 
because the nuclear distributions do not vanish at x = 1 . 
However, it is more practical at this stage to use the 
form due to the lack of large-x data. In any case, the 
large x (> 0.8) region does not contribute to the anomaly 
because of the NuTeV kinematics. For investigating the 
difference between the modifications, we define a function 



£v(x) 



w di] (x, Q 2 , A, Z) - w u „ (x, Q 2 , A, Z) 
w dv (x, Q 2 , A, Z) + w Uv (x, Q 2 , A, Z) 



(3) 



The function e v depends also on A, Z, and Q 2 , but these 
factors are abbreviated for simplicity. 

Writing the cross sections in terms of the nuclear par- 
ton distribution functions (PDFs) in the leading order of 
a s , we obtain a modified PW relation with small correc- 
tion factors. Expanding the modified relation in terms 



of the correction factors and retaining only the leading 
correction of e v , we obtain 
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+ 0(e 2 v ) + 0(e n ) + 0(e s ) + 0(e c ). 



(4) 



Here, 0{e) indicates a correction of the order of e, and 
the variable y is given by the energy transfer q° and 
neutrino energy E u as y = q°/E v . The derivation of 
Eq. Q is found in Ref. [TJ. The details are investi- 
gated for the first correction factor e v (x) in this paper. 
Other corrections come from neutron excess, strange- 
antistrange asymmetry (s — s), and charm-anticharm 
asymmetry (c — c). They are denoted 0(e n ), 0(e s ), 
and 0(e c ) in Eq. (jJJ. The e factors are defined by 



(N-Z){u v -d v )/[A{u v +d v )], 
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= c^/[w v (u v + d v )] with w v = (w<i v +w Uv )/2 and 
q A — q A . The correction factor e v is related to 
the isospin violation as shown in Sec. IIV Bl The studies 
of these contributions to sin 2 6w should be found else- 
where. In deriving Eq. @, differential cross sections are 
used instead of the total ones for taking the PW ratio. 
The equation is useful for finding theoretical possibilities; 
however, it is practically limited for estimating a numer- 
ical effect because NuTeV kinematical conditions should 
be taken into account. This point is discussed in Sec. 
IIVI Investigating the modified PW relation, we find that 
(2) the e v (x) factor contributes to the sin 2 6\y measurement 



in the neutrino scattering. We explain the distribution 
e v (x) and its effect on the sin 2 By/ determination in the 
following sections. 



III. x ANALYSIS METHOD FOR 
DETERMINING e v (x) 

The nuclear modification difference between the dis- 
tributions u v and d v is not known at this stage. There 
are few theoretical guidelines for calculating the differ- 
ence. Of course, the valance-quark distributions are con- 
strained by the baryon-number and charge conservations, 
so that there should be some restrictions to their nu- 
clear modifications. From this consideration, two possi- 
bilities were proposed for e v (x) in Ref. However, 
the conservation conditions are not enough to impose in- 
dependent shape of these modifications. It is, therefore, 
an appropriate way is to show the modification difference 
by analyzing available experimental data. 

We express the difference by a number of parameters, 
which are then determined by a x 2 analysis of nuclear 
data. The difference between the valence-quark modifi- 
cations is denoted 

Aw v (x, Q 2 ,A, Z) = w Uv (x,Q 2 ,A,Z) - w dv (x,Q 2 , A, Z), 

(5) 
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and their average is w v as used in the previous section. 
The difference at Q 2 —l GeV 2 (= Q 2 ,) is expressed by 
four parameters, a' v , b' v , d v , and d' v : 



Aw v {x,Ql,A,Z) 
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a' v (A, Z) + b' v x + c' v x 2 + d' v x 3 
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For the functions w v , Wq, and ui g , which indicate nuclear 
modifications of average valence-quark, antiquark, and 
gluon distributions, we use the ones obtained by a re- 
cent global analysis in Ref. J19j . The reason of selec ting 
the functional form of Eq. © is discussed in Ref. |20| . 
We briefly explain the essential point. The A depen- 
dence of the modification is assumed to be proportional 
to 1 — 1/A 1 / 3 by considering nuclear volume and surface 
contributions to the cross section [2l| . It should be noted 
that such A dependence could be too simple to describe 
the distributions. The factor 1/(1 — x)' 3 " is introduced 
so as to explain the large- a; Fermi- motion part. Looking 
at typical data of F^jF® ratios, we find that a cubic 
function seems to be appropriate for the x dependence. 
However, we should aware that an appropriate functional 
form is scarcely known for Aw v . 

The parameters a' v , b' v , c' v , and d' v are determined by 
analyzing experimental data for structure-function ratios 
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and Drell-Yan cross-section ratios a v DY la 



_pA I pA' 



The details of the experimental data sets are discussed 
in Ref. 0]. The F£ /F^' data are from the European 
Muon Collaboration (EMC) the SLAC-E49, E87, 
E139, and E140 Collaborations^! . the Bologna-CERN- 
Dubna-Munich-Saclay (BCDMS) Collaboration 0, the 
New Muon Collaboration (NMC) ;25|, the Fermilab-E665 
Collaboration [H, and HERMES Collaboration [27|. 
The Drell-Yan data are taken by the Fermilab-E772 and 
E866/NuSea Collaborations |28j. These data are for the 
nuclei: deuteron, helium-4, lithium, beryllium, carbon, 
nitrogen, aluminum, calcium, iron, copper, krypton, sil- 
ver, tin, xenon, tungsten, gold, and lead. The total num- 
ber of the data is 951. Because the small- a; data are taken 
in the small Q 2 region, our leading-twist analysis could 
be affected especially at small x |19| . 

One of these parameters is fixed by baryon-number or 
charge conservation, and a' v is selected for this fixed quan- 
tity. In addition, the parameter f3 v is taken as the same 
value in the analysis of Ref. 0] ((3 V = 0.1). Therefore, 
there are three free parameters in the analysis. The total 
X 2 is defined by 
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(7) 



where Rj indicates the ratios, 



F 2 A /Ff and a£i/<y 



The theoretical ratios i?* are calculated by evolving 
the nuclear PDFs with the valence-quark modification 
difference in Eq. © by the Dokshitzcr-Gribov-Lipatov- 
Altarelli-Parisi (DGLAP) equations to experimental Q 2 



points. The total x 2 is minimized by the subroutine 
MINUIT [29j. Running the subroutine, we obtain opti- 
mized values for the parameters and also a Hessian ma- 
trix for error estimation. The uncertainty of the weight 
function Aw v (x, Q 2 , £) is calculated with the Hessian ma- 
trix: 

L?A I r>2M2 A - \ " ( dAw v (x,Q 2 ,£,) 

[SAw v (x,Q )} =Ax }^[ gt. 
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where a parameter is denoted and the optimized point 
is denoted £. The A\ 2 value is taken so that the confi- 
dence level is the one-cr-error range for the normal dis- 
tribution |19|. Because the parameter number is three, 
we have A\ 2 = 3.527. This uncertainty estimation is es- 
pecially important in discussing an effect on the NuTeV 
anomaly. 



IV. RESULTS 

A. Valence-quark modification difference 
between u v and d v 

We explain x 2 analysis results. From the analysis, we 
obtain Xmin/ d-°-f -=L57. The parameters are defined at 
Q 2 =l GeV 2 , and their optimum values obtained by the 
fit are shown in Tabled We note that the constant a' v 
is fixed by one of the baryon-number and charge con- 
servations, and it depends on nuclear species. Several 
nuclei (D, C, Ca, Fe, Ag, Au) are selected for listing 
its values. It is almost constant from a small nucleus 
K(-D) = 0.00612) to a large one {a' v (Au) = 0.00615). 
The parameters b' v , c' v , and d' v are determined by the 
analysis. We notice huge errors which are almost an or- 
der of magnitude larger than the optimum parameter val- 
ues. This fact indicates that the determination of Aw v is 
almost impossible at this stage. However, it is important 
to show an uncertainty range of Aw v in order to compare 
with the NuTeV anomaly. 



TABLE I: Parameters obtained by the analysis. The constant 
a' v depends on nuclear species. Its value is listed for typical 
nuclei. 



parameter 


value 


b' v 


-0.0905 ± 0.4044 


c' 


0.239 ± 1.300 


d' v 


-0.193 ± 1.126 


a v {D) 
a' v (C) 
a' v (Ca) 
a' v {Fe) 
<{Ag) 
a' v (Au) 


0.00612 
0.00612 
0.00612 
0.00613 
0.00614 
0.00615 
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FIG. 1: (Color online) Nuclear modification difference Aw v = 
W Uv — Wd v between u v and d v for the iron nucleus. The curve 
indicates the obtained Aw v at Q 2 =l GeV 2 , and the shaded 
area shows the uncertainty range calculated by the Hessian 
method. 



The difference w Uv — Wd v is plotted as a function of x at 
Q 2 = l GeV 2 in Fig. Q]by using the optimum parameters. 
The solid curve is the distribution obtained by the % 2 
analysis. The shaded area is the one-tr uncertainty range 
which is calculated by the Hessian method in Eq. JSJ. 
Although the distribution is positive at small x and it 
becomes negative at large x, its functional form is not 
clear if the uncertainties are considered. Because the 
uncertainty band is much larger than the distribution 
itself, we obviously need future experimental efforts to 
find an accurate distribution. 



B. Effect on sin 2 6 W 

The average Q 2 of the NuTeV neutrino and antineu- 
trino scattering experiments is about Q 2 =20 GeV 2 . The 
valence-quark distributions in the iron nucleus are ob- 
tained at Q 2 =l GeV 2 in the previous section. The nu- 
clear PDFs as well as the nucleonic PDFs are evolved 
to the ones at Q 2 =20 GeV 2 by the DGLAP evolution 
equations. Then, the nuclear modifications w Uv and 
are calculated at Q 2 =20 GeV 2 by using Eq. J2J, and the 
obtained function £^(2;) is shown in Fig. [5] 

The solid curve is the distribution obtained by the 
current analysis, and its uncertainties are shown by the 
shaded area. The previous results in Ref. are a l so 
shown in the figure. These distributions are proposed 
as the ones which satisfy the baryon-number and charge 
conservations. The dashed curve indicates the distribu- 
tion by the model-1, in which the conservation integrands 
are assumed to vanish in the leading order of the correc- 
tion factors. The dotted curve indicates the distribution 
by the model-2, in which the global \ 2 analysis results for 
the nuclear PDFs in Ref. [2(j are used. Three curves are 
much different; however, they are certainly within the 
error band. It suggests that they should be consistent 
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FIG. 2: (Color online) The function e v (x) at Q 2 =20 GeV 2 
for the iron nucleus. The solid curve indicates the current 
analysis result, and the shaded area shows the uncertainty 
range. The dashed and dotted curves indicate the model-1 
and model-2 results, respectively, in Ref. |ll| . 



within the uncertainties. 

Because a Monte-Carlo code, instead of the PW rela- 
tion, is used in the experimental analysis for extracting 
sin 2 9w, it is not straightforward to compare theoretical 
results with the NuTeV deviation. Obviously, NuTeV 
kinematical conditions should be taken into account. In 
particular, few data exist in the large- x region, where 
e v (x) is large. Fortunately, such kinematical factors are 
provided in Ref. [T^- In order to use them, our PDFs 
should be related to the NuTeV convention |3(j • The up- 
and down- valence quark distributions per nucleon are de- 
fined by 



Zxu v + Nxd v Zu 



Nv* 

vp ' vn 



A 



A 



A Zxd v + Nxu v Zd* + Nd* v 

xd v = w dv - 



A 



A 



(9) 



where the middle expressions with w Uv and Wd v are our 
notations, and the PDFs with the asterisk (*) are the 
NuTeV expressions. Using Eq. @, we find that the 
distribution e v (x) is related to the isospin- violating dis- 
tributions: 



5u* v 
5d* 



d* vn = —e v (w Uv + w d Jx u v , 

u* vn = +e v (w Uv + Wd v )x d v . (10) 



Therefore, according to the NuTeV convention of the nu- 
clear PDFs, we have been investigating isospin-violation 
effects in the nucleon and their nuclear modifications by 
the function £^(3;). However, we should aware that such 
isospin-violation effects could possibly include other nu- 
clear effects which may not be related to the isospin vio- 
lation. In any case, using the NuTeV functionals in Ref. 



4 



[13|. we need to calculate the integral 



A (sin 2 6 W ) 



dx { F[Sul, x] 5u*(x) 
+ F[Sd* v ,x] 5d* v (x)}, 



for estimating a correction to the NuTeV measurement. 
All the NuTeV kinematical effects are included in the 
functionals F \6u *,, x ] and F[Sd*, x] which are provided in 
Fig. 1 of Ref. [13. It should be noted that our sign 
convention of the correction A (sin 2 Qy/) is opposite to 
the NuTeV one. 

Using Eqs. (|10|l and together with the distribu- 
tion e v (x) obtained in Sec. IIV Al we calculate the effect 
on the NuTeV value. The results are shown in Fig. 
as a function of Q 2 . The deviation A(sin 2 9w) is shown 
by the solid curve, and the shaded area corresponds to 
the one-er-error range. For comparison, the NuTeV de- 
viation 0.0050 is shown by the dashed line in the figure. 
We notice that the correction is not strongly dependent 
on Q 2 . Even if the uncertainty range is considered, the 
contribution is smaller than 0.0050. In fact, calculating 
the correction at Q 2 =20 GeV 2 , which is approximately 
the average Q 2 of NuTeV measurements, we obtain 



A(sin 2 6 W ) = 0.0004 ± 0.0015. 



that the uncertainty estimation depends on the analysis 
conditions. For example, a certain A and x dependent 
functional form is assumed in the y 2 analysis and such 
uncertainties are not included in Fig. Eland Eq. Ijl2|l . 
The situation is the same as the nucleonic PDF case, for 
example in Ref. |3l| , where such functional uncertainties 
are not also estimated. In addition, the uncertainties due 
to other parameter errors in Ref. 0] are not included. 
There is a possibility that the uncertainty 0.0015 could 
be underestimated. 

Because of these issues, it is, strictly speaking, too 
early to exclude the mechanism for explaining the NuTeV 
anomaly. However, if the NuTeV anomaly were to be ex- 
plained solely by the nuclear modification of the PDFs, 
the magnitude should be a factor of ten larger than the 
global analysis result. Furthermore, the obtained error 
is a factor of three smaller than the NuTeV discrepancy. 
Therefore, it is unlikely at this stage that the anomaly 
is explained by such an effect according to the analy- 
sis of nuclear data. In any case, because the difference 
W Uv — u>d v cannot be determined from current experi- 
mental data, we need future experimental efforts to find 
it. Because it is related to nuclear valence-quark distri- 
butions, possibilities are neutrino scattering experiments 
at future neutrino facilities such as MINERVA 32] and 
(12) neutrino factories 1331 . 
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V. SUMMARY 

We have extracted the difference between nuclear mod- 
ifications of u v and d v by the x 2 analysis of nuclear F2 
and Drell-Yan data. We found that it is difficult to de- 
termine it at this stage due to the lack of data which 
are sensitive to the difference. Such a difference con- 
tributes to the NuTeV determination of sin 2 9w because 
the iron target was used in the experiment. Taking the 
NuTeV kinematics into account, we estimated the effect 
on the sin 2 6y/ . It is not large enough to explain the 
whole NuTeV deviation. However, the nuclear modifica- 
tion difference cannot be accurately determined at this 
stage, we need further efforts to find it and its effect on 
the NuTeV determination of sin 2 9w ■ 



FIG. 3: (Color online) Effects on the NuTeV sin 2 9 W deter- 
mination as a function of Q 2 by taking into account the ex- 
perimental kinematics. The curve indicates the \ 2 analysis 
result, and the shaded area shows the uncertainty range. The 
NuTeV deviation 0.0050 is shown by the dashed line. 



It indicates that the whole NuTeV anomaly could 
not be explained by the nuclear modification difference 
W Uv — Wd v at this stage. However, we should be careful 
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